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Abstract. We study the exterior derivative as a symmetric unbounded operator on square 
integrable 1-forms on a 3D bounded domain D. We aim to identify boundary conditions that render 
this operator self-adjoint. By the symplectic version of the Glazman-Krein-Naimark theorem this 
amounts to identifying complete Lagrangian subspaces of the trace space of fJ"{curl, D) equipped 
with a symplectic pairing arising from the A-product of 1-forms on dD. Substantially generalizing 
earlier results, we characterize Lagrangian subspaces associated with closed and co-closed traces. In 
the case of non-trivial topology of the domain, different contributions from co-homology spaces also 
distinguish different self-adjoint extension. Finally, all self-adjoint extensions discussed in the paper 
are shown to possess a discrete point spectrum, and their relationship with curl curl-operators is 
discussed. 
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1. Introduction. The curl operator is pervasive in field models, in particular 
in electromagnetics, but hardly ever occurs in isolation. Most often we encounter a 
curl curl operator and its properties are starkly different from those of the curl alone. 
We devote the final section of this article to investigation of their relationship. 

The notable exception, starring a sovereign curl, is the question of stable force- 
free magnetic fields in plasma physics. They are solutions of the eigenvalue problem 

aeM\{0}: curlH = aH, (1.1) 

posed on a suitable domain, see [11,20,26,33]. A solution theory for (|1.1|) must scru- 
tinize the spectral properties of the curl operator. The mature theory of unbounded 
operators in Hilbert spaces is a powerful tool. This approach was pioneered by R. 
Picard [31,34,35], see also [39]. 

The main thrust of research was to convert curl into a self-adjoint operator by 
a suitable choice of domains of definition. This is suggested by the following Green's 
formula for the curl operator: 

/ curlu • V — curlv ■ udx = / (uxv) ndS', (1-2) 

Jd JdD 

for any domain D C M.^ with sufhciently regular boundary dD and u, v G C^{D). 
This reveals that the curl operator is truly symmetric, for instance, when acting on 
vector fields with vanishing tangential components on dD. 

On bounded domains D several instances of what qualifies as a self-adjoint curl 
operators were found. Invariably, their domains were defined through judiciously cho- 
sen boundary conditions. It also became clear that the topological properties of D 
have to be taken into account carefully, see [34, Thm. 2.4] and [39, Sect. 4]. 

In this paper we carry these developments further with quite a few novel twists: 
we try to give a rather systematic treatment of different options to obtain self-adjoint 
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curl operators. It is known that the curl operator is an incarnation of the exterior 
derivative of 1-forms. Thus, to elucidate structure, we wiU mainly adopt the perspec- 
tive of differential forms. 

Further, we base our considerations on recent discoveries linking symplectic alge- 
bra and self-adjoint extensions of symmetric operators, see [16] for a survey. In the 
context of ordinary differential equations, this connection was intensively studied by 
Markus and Everitt during the past few years [14]. They also extended their investiga- 
tions to partial differential operators like A [15]. We are going to apply these powerful 
tools to the special case of curl operators. Here, the crucial symplectic space is a 
Hilbert space of 1-forms on dD equipped with the pairing 



We find out, that it is the Hodge decomposition of the trace space for 1-forms on 
D that allows a classification of self-adjoint extensions of curl: the main distinction 
is between boundary conditions that impose closed and co-closed traces Moreover, 
further constraints are necessary in the form of vanishing circulation along certain 
fundamental cycles of dD. This emerges from an analysis of the space of harmonic 
1-forms on dD as a finite-dimensional symplectic space. For all these self-adjoint curls 
we show that they possess a complete orthonormal system of eigenfunctions. 

In detail, the outline of the article is as follows: The next section reviews the con- 
nection between vector analysis and differential forms in 3D and 2D. Then, in the third 
section, we introduce basic concepts of symplectic algebra. Then we summarize how 
those can be used to characterize self-adjoint extensions through complete Lagrangian 
subspaces of certain factor spaces. The fourth section applies these abstract results to 
trace spaces for 1-forms and the corresponding exterior derivative, that is, the curl 
operator. The following section describes important complete Lagrangian subspaces 
spawned by the Hodge decomposition of 1-forms on surfaces. The role of co-homology 
spaces comes under scrutiny. In the sixth section we elaborate concrete boundary con- 
ditions for self-adjoint curl operators induced by the complete Lagrangian subspaces 
discussed before. The two final sections examine the spectral properties of the classes 
of self-adjoint curls examined before and explore their relationships with curl curl 
operators. Frequently used notations are listed in an appendix. 

2. The curl operator and differential forms. In classical vector analysis 
the operator curl is introduced as first order partial differential operator acting on 
vector fields with three components. Thus, given a domain Z3 C M'^ we may formally 
consider curl : Cjf (Z3) ^ C^{D) as an unbounded operator on L^{D). Integration 
by parts according to (|1.2[) shows that this basic curl operator is symmetric, hence 
closable [38, Ch. 5]. Its closure is given by the minimal curl operator 



Note, that curlmax is no longer symmetric, and neither operator is self-adjoint. This 
motivates the search for self-adjoint extensions curlg : I?(curls) C L'^{D) i-^ L^{D) 
of curlijiin. If they exist, they will satisfy, c.f. [16, Example 1.13], 




cxiyU,^ : Ho{cnv\,D) ^ L\D) . 
Its adjoint is the maximal curl operator, see [34, Sect. 0], 
curlniax curl^in : if (curl, D) ^ 



(D). 



(2.1) 



(2.2) 



curl] 



■min 



C curls C curl 



■max 



(2.3) 
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Remark 1. The classical route in the study of self-adjoint extensions of symmetric 
operators is via the famous Stone- von Neumann extension theory, see [38, Ch. 6]. It 
suggests that, after complexification, we examine the deficiency spaces 



Lemma 2.1. The deficiency spaces from p.4p satisfy dim — oo. 

Proof. Let : R'^ \ {0} i-^ C'^'^ be a fundamental solution (dyad) of curl ±i, that 
is, G = (curl — V^V)I$, where $(a;) — exp(-|a;|)/(47r|a;|) is the fundamental 

solution of -A + 1, and V (J^, Then, for any ip e C-(R3)|g^, 



satisfies u G i?(curl, D) and curlu ±m~Q.U 

From Lemma l^Tj we learn that reveal little about the structure governing self- 
adjoint extensions of curl. Yet the relationship of curl and differential forms suggests 
that there is rich structure underlying self-adjoint extensions of curl,„in. 

2.1. Differential forms. The curl operator owes its significance to its close 
link with the exterior derivative operator in the calculus of differential forms. We 
briefly recall its basic notions and denote by M an m-dimensional compact orientable 
manifold with boundary dM . If Af is of class it can be endowed with a space of 
differential forms of degree k, < k < m: 

Definition 2.2 (Differential fc-form). A differential form of degree k (in short, a 
k-form) and class C\ I G Nq, is a C'^ -mapping assigning to each x ^ M an alternating 
k-multilinear form on the tangent space Tx{M). We write A*''''(M) for the vector space 
of k-forms of class C' on M, A'''\M) — {0} for k <0 or k > m. 

Below, we will usually drop the smoothness index I, tacitly assuming that the 
forms are "sufficiently smooth" to allow the respective operations. 

The alternating exterior product of multilinear forms gives rise to the exterior 
product A : A'^(M) x A-'(Af) ^ A'''+^(Af) by pointwise definition. We note the graded 
commutativity rule u! Arj — {—l)^^rif\u) for a; £ A*^, 7] G A^ . Further, on any piecewise 
smooth orientable fc-dimensional sub-manifold of M we can evaluate the integral J^^ u 
of a fc-form lo over a fc-dimensional sub-manifold E of M [10, Sect. 4]. 

This connects to the integral view of fc-forms as entities that describe additive 
and continuous (w.r.t. to a suitable deformation topology) mappings from orientable 
sub-manifolds of M into the real numbers. This generalized differential forms are 
sometimes called currents and are studied in geometric integration theory [13,17]. 
From this point of view differential forms also make sense for non-smooth manifolds. 

From the integral perspective the transformation (puUback) ^*lo of a fc-form 
under a sufficiently smooth mapping <^ : M ^ M appears natural: $*aj is a fc-form 
on M that fulfills 



A/'(cUrlniax ±« • Id) C I?(curlmax) • 



(2.4) 





(2.5) 



for all fc-dimensional orientable sub-manifolds E of M. We remark that puUbacks 
commute with the exterior product. 
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If i : dM i-^ M stands for the inclusion map, then the natural trace of a fc-form 
Lo G A'=(M) on dM is defined as i*Lu. It satisfies the following commutation relations 

i*{u; Arj) ^ {i*uj) A{i*r]) and d {i* uj) ^ i* {d uj) . (2.6) 

The key operation on differential form is the exterior derivative 

d : A'=(M) ^ A'=+i(M) , (2.7) 

which is connected with integration through the generalized Stokes theorem 

( duj= ( uj yuieA'^'^iM) (2.8) 

and all orientable piecewise smooth sub-manifolds of M. In fact, (|2.8p can be used to 
define the exterior derivative in the context of geometric integration theory. This has 
the benefit of dispensing with any smoothness requirement stipulated by the classical 
definition of d. We have = and, obviously, (|2.8p and (|2.7p imply $* o d = d o$*. 
Since one has the graded Leibnitz formula 

d{ujAr]) ^ d Lu A +{-!)''' Lu Ad ri e A'=(A.f), 77 G A^(Af) (2.9) 

exterior derivative and exterior product enter the crucial integration by parts formula 

I dujAr]+{-l)'' [ ujAd7]= [ i*ujAi*T] Vuj E A'' {M) , rj e (M) , (2.10) 

where S is an orientable sub-manifold of M with dimension k + j + 1 and canonical 
inclusion i : S dT,. 

2.2. Metric concepts. A metric g defined on the manifold M permits us to 
introduce the Hodge operator -kg : A^{M) \-* A"^~^{M). It gives rise to the inner 
product on A^{M) 

{uj,r])k,M / t^A*gr/, uj,f] e A'^iM) . (2.11) 

Thus, we obtain an L^-type norm ||-|| on A'^(M). Completion of smooth fc-forms with 
respect to this norm yields the Hilbert space L^(A'^(M)) of square integrable (w.r.t. 
g) fc-forms on M. Its elements are equivalence classes of fc-forms defined almost every- 
where on M. Since Lipschitz manifolds possess a tangent space almost everywhere, 
for them L^(A'^(M)) remains meaningful. As straightforward is the introduction of 
"Sobolev spaces" of differential forms, see [1, Sect. 1], 

W''{d,M) {uj e L^{A^{M)) : duj E L^{A^{M))} , (2.12) 

which are Hilbert spaces with the graph norm. The completion of the subset of smooth 
fc-forms with compact support in W''{d,M) is denoted by IVf (d, Af). 
By construction, the Hodge star operator satisfies 

= (-l)(™-fc)fc Id . (2.13) 

Now, let us assume dM = 0. Based on the inner product (|2.1ip we can introduce the 
adjoint d* : VF'^+^(d,M) W''{d,M) of the exterior derivative operator by 

{dLo,7j)k+i,M ^ {Lu,d*rj)k,M VLueW''+\d,M),r^eWS{d,M) , (2.14) 
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and an explicit calculation shows that 

d* = ^ d * : A'^'+i ^ A'^' . (2.15) 

Furthermore, d^ = implies (d*)^ = 0. Eventually, the Laplace-Beltrami operator is 
defined as 

Am = dd*+d*d:A'^' — > A'^'. (2.16) 

2.3. Vector proxies. Let us zero in on the three-dimensional "manifold" D. 
Choosing bases for the spaces of alternating /c-multilinear forms, differential fc-forms 
can be identified with vector fields with (^) components, their so-called "vector 
proxies" [1, Sect. 1]. The usual association of "Euclidean vector proxies" in three- 
dimensional space is summarized in Table [^?T] The terminology honours the fact that 
the Hodge operators ★ : K^{D) i-» h^{D) and A°(D) i-^ K^{D) connected with the 
Euclidean metric of 3-space leave the vector proxies invariant (this is not true in 2D 
since *^ = — 1 on 1-forms). In addition the exterior product of forms is converted into 
the pointwise Euclidean inner product of vector fields. Thus, the inner product (•, ■)k,D 
of fc-forms on D becomes the conventional L^(D) inner product of the vector proxies. 
Further, the spaces W^{A^D) boil down to the standard Sobolev spaces H^{D) (for 
fc ^ 0), ff(curl, D) (for fc = 1), i?(div, D) (for k = 2), and L'^{D) (for fc = 3). 



Differential form uj Related function u or vector field u 

X ^ oj{x) u{x) := uj{x) 

X {v 1-^ uj{x)(\-)} u(a;) ■ v := aj(a;)(v) 

X {(vi, V2) ^ a)(a;)(vi, V2)} u{x) ■ (vi x V2) a'(a;)(vi, V2) 

X ^ {(vi,V2,V3) w(a;)(vi, V2, V3)} ii(x) det(vi, V2, V3) := w(a;)(vi, V2, V3) 

Table 2.1 

The standard choice of vector proxy u,u for a differential form lo in M.^ . Here, ■ denotes the 
Euclidean inner product of vectors inM.^, whereas X designates the cross product. 



Using Euclidean vector proxies, the curl operator turns out to be an incarnation 
of the exterior derivative for 1-forms. More generally, the key first order differential 
operators of vector analysis arise from exterior derivative operators, see Figure 12.11 
Please note that, since the Hodge operator is invisible on the vector proxy side, curl 
can as well stand for the operator 

curl < — > *d : Ai(D) Ai(D) , (2.17) 

which is naturally viewed as an unbounded operator on L^{A^{D)). Thus, (|2.17|) puts 
the formal curl operator introduced above in the framework of differential forms on 
D. 

Translated into the language of differential forms, the Green's formula p.2p be- 
comes a special version of (|2.10p for fc = j = 1. However, due to (|2.13p . p.2p can also 
be stated as 

i*dLJ,7j)i,D~icJ,*dr])i,D^ [ i*iJAi*v, Lu,veW\d,D). (2.18) 

JdD 

A metric on induces a metric on the embedded 2-dimensional manifold dD. 
Thus, the Euclidean inner product on local tangent spaces becomes a meaningful 
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0-forms A"(i:>) 



1-forms A^{D) 



* d * 



functions D ^ R 
grad grad* = div 

field intensity vector fields D t-^ M.^ 



d* =*d* 



curl 



curl* = curl 



2-forms A^{D) 



flux density vector fields D 



★ d* 



div 



div* 



grad 



3-forms A^{D) scalar densities ^ R 

Fig. 2.1. Differential operators and their relationship with exterior derivatives 



Differential forms 


Related function u or vector field u 


X 1-^ uj(x) 


u(x) :— uj{x) 


X 1-^ {v Lli{x){\-)} 


u{x) ■ V := U!{x){\-) 


X ^ {(vi,V2) ^ tj(a;)(vi, V2)} 


u{x) det(vi, V2, n(a;)) := w(a;)(vi, V2) 



Table 2.2 

Euclidean vector proxies for differential forms on dD. Note that the test vectors v, vi , V2 have 
to be chosen from the tangent space T^idD). 



concept and Euclidean vector proxies for fc-forms on dD, fc = 0, f , 2, can be defined 
as in Table O see Table O 

This choice of vector proxies leads to convenient vector analytic expressions for 
the trace operator i*: 



u e A"{D) 

Lo e A^D) 
UJ e A^D) 



I UJ 

i*uj 
i*uj 
i*uj 



"fu{x) 
-ftu{x) 

0, 



— u{x), u : D I 

— u{x) — {u{x) ■ n(a;))n(a;), u : D 
= u{x) ■ n{x), u : D 



where x G dD. Further, the customary vector analytic surface differential operators 
realize the exterior derivative for vector proxies, see Figure [2?2l 

Remark 2. Vector proxies offer an isomorphic model for the calculus of differential 
forms. However, one must be aware that the choice of bases and, therefore, the de- 
scription of a differential form by a vector proxy, is somewhat arbitrary. In particular, 
a change of metric of space suggests a different choice of vector proxies for which 
the Hodge operators reduce to the identity. Thus, metric and topological aspects are 
hard to disentangle from a vector analysis point of view. This made us prefer the 
differential forms point of view in the remainder of the article. 



3. Self-adjoint extensions and Lagrangian subspaces. First, we would like 
to recall some definitions of symplectic geometry. Then, we will build a symplectic 
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O-forms A°{dD) 



functions 



d* = — * d * 



l-forms A^{dD) 



gradg 



grad* = —diva 



vector fields 



d* = — * d * 



2-fornis A'^{dD) 



curia 



curlg = curia = n x gradg 



scalar densities 



Fig. 2.2. Exterior derivative for Euclidean vector proxies on 2-manifolds 



space associated to a closed symmetric operator. The reader can refer to [15, 16] for a 
more detailed treatment. 

3.1. Concepts from symplectic geometry. Symplectic geometry offers an 

abstract framework to deal with self-adjoint extensions of symmetric operators in 
Hilbert spaces. Here we briefly review some results. More information is available 
from [28]. 

Definition 3.1 (Symplectic space). A real symplectic space S is a real linear 
space equipped with a symplectic pairing [■,■] (symplectic bilinear form, symplectic 
product) 



[■,■]: SxS — 
[ami + a2U2,v] 
[u,v] = -[v,u], 
[ [u,S] = =^ 



ai wi,?; + a2\u2,v\ 



u = 



(linearity) 
(skew symmetry) 
( non- degeneracy ) 



(3.1) 



Definition 3.2. Let L be a linear subspace of the symplectic space S 
(i) The symplectic orthogonal of L is L'^ = {u € S : [u, L] = 0}; 
(ii) L is a Lagrangian subspace, if L C i. e. [u,v] = for all u and v in L; 
(Hi) A Lagrangian subspace L is complete, if = L. 

In the case of finite dimensional symplectic spaces, symplectic bases offer a con- 
venient way to build complete Lagrangian subspaces, see [14, Example 2]. 

Definition 3.3. Let {S, [■,■]) be a real sym.plectic space with dimension 2n (the 
dimension has to be even so that the pairing [■, ■] can be non degenerate). A symplectic 
basis of S is a basis of S satisfying 



[Ui,Uj\ 



Jij with J = 



Onxn Inxn 
~Inxn O^xn 



(3.2) 



Simple linear algebra proves the existence of such bases: 

Lemma 3.4. For any symplectic space with finite dimension 2n, there exists a 
(non unique) symplectic basis. 

Remark 3. As soon as we have found a symplectic basis {uj}?^!, it provides many 
complete Lagrangian subspaces 

• the n first vectors of a symplectic basis span a complete Lagrangian 

subspace. 
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• the n last vectors {ui}f2n+i '^^ ^ symplectic basis span a complete Lagrangian 
subspace. 

• for any a : |l,n] ^ |0, 1], {ui+CT(i)n}"=i is a complete Lagrangian subspace. 
We recall some more facts about finite dimensional symplectic spaces 
Lemma 3.5. Every complete Lagrangian subspace of a finite dimensional sym- 
plectic space S of dimension 2n is n- dimensional. Moreover, it possesses a basis that 
can be extended to a symplectic basis of S . 

3.2. Application to self-adjoint extensions of a symmetric operator. 

Let iJ be a real Hilbert space and T a closed symmetric linear operator with dense 
domain X'(T) C H . We denote by T* its adjoint. Let us first recall, see [38], that each 
self-adjoint extension of T is a restriction of T*, which is classically written as 

T c C T*. (3.3) 

Next, introduce a degenerate symplectic pairing on X'(T*) by 

[•,•]: 2?(T*) X P(T*) — >R such that [u,v] = (T*u,v) - (u,T*v). (3.4) 

From the definition of T*, the symmetry of T, and the fact T** = T, we infer that, 
see [14, Appendix], 

[u + uq,v + vo\ = [u,v], V-uo,'yo e ©(T), Vu, u e P(T*), 
M e P(T*), [w,w] = 0, VwgP(T*), =^ weP(T). 

As a consequence, we obtain a symplectic factor space, see Appendix of [14], 

Lemma 3.6. The space S = ^X'(T*)/X'(T), [•, ■]j is a symplectic space. The graph 

norm on P(T*) induces a factor norm on S and, due to (|3.5p . the symplectic pairing 
[•, •] is continuous with respect to this norm 

|[u,v]|2 < {\\uf + \\J*uf) ■ {\\vf + \\J*vf) Vu G V(J*), V e V{T*) , 

Let L ® T>{T) denotes the pre-image of L under the factor map I?(T*) i-^ S. 

Corollary 3.7. The symplectic orthogonal V'^ of any subspace V of S is closed 
(in the factor space topology). Any linear subspace L oi S defines an extension of 
T through 

TcTl :=T*|ieP(T) CT* . (3.6) 

This relationship allows to characterize self-adjoint extensions of T by means of the 
symplectic properties of the associated subspaces L. This statement is made precise 
in the Glazman-Krein-Naimark Theorem, see Theorem 1 of [14, Appendix]. 

Theorem 3.8 (Glazman-Krein-Naimark Theorem symplectic version). The map- 
ping L Ti is a bisection between the space of complete Lagrangian subspaces of S 
and the space of self- adjoint extensions ofT. The inverse mapping is given by 

L = V{Jl)/V{T) . (3.7) 
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4. Symplectic space for curl. Evidently, the unbounded curl operators intro- 
duced in Section [5] (resorting to the vector proxy point of view) fits the framework of 
the preceding section and Theorem l3.8l can be apphed. To begin with, from (|2.ip and 
(|2.2p we arrive at the symplectic space 



^curi //(curl, D)/Ho{curl D) . (4.1) 
By (jl.2p it can be equipped with a symplectic pairing that can formally be written as 



\aD 



f (u(y) xv(y)).n(y)d5(y) , (4.2) 

JdD 



for any representatives of the equivalence classes of S'curi- From (j4.ip it is immediate 
that 5* is algebraically and topologically isomorphic to the natural trace space of 
//(curl,!?). 

By now this trace space is well understood, see the seminal work of Paquet [30] 
and [6-9] for the extension to generic Lipschitz domains. To begin with, the topology 
of S'curi is intrinsic, that is, with D' := M!^ \ D, the norm of 

S^^,i := //(curl, i?')/-ffo(curl, D') (4.3) 

is equivalent to that of S'curi! both spaces are isomorphic algebraically and topologi- 
cally. This can be proved appealing to an extension theorem for //(curl, D). Moreover, 
the pairing [•, -J^^ identifies 5curi with its dual S'curi- 

Lemma 4.1. The mapping Scuri > S'curi; u i— > {v i-^ [u, v]g^} is an isomorphism. 

Proof. Given u G Scuri, let u £ //(curl, D) solve 

curl curl u + u = in D , 7tu = u on 51? . (4-4) 
Set V := curlu £ //(curl, D) and v := -jtv e Scuri- By pT^ 



u,v]gjj= / curl u • V — curl V • udx 



ID 

^|cUrlu|2-H|u|2dx= ||v||jj(^„rl,D) l|u|lff(curl,D) > IM I 5<,„, J I ^ II S,„,, , 

as I|v||ff(curi,i5) = ll"ll/f(curi,L>)- Wc immediately conclude 

veSe„., l|V|ISe„, 

□ 

The trace space also allows a characterization via surface differential operators. It 

relies on the space H^^ (dD) of tangential surface traces of vector fields in {H^ (D))^ and 
— i 1 

its dual H( ^ (dD) :— (H.^ (dD))'. Then one finds that, algebraically and topologically, 
Scuri is isomorphic to 

S'curi = H"5(curla,aZ?) -.^ {v e U;^^D) : curia v G H-^dD)} . (4.5) 

The intricate details and the proper definition of curl^ can be found in [9] . 

When we adopt the perspective of differential forms, the domain of curlmax is 
the Sobolev space H^{d, D) of 1-forms. Thus, S'curi has to be viewed as a trace space 
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of 1-forms, that is, a space of 1-forms (more precisely, 1-currents) on dD. In analogy 
to (|2.12p and (|4.5p we adopt the notation 

S,^ri^W-"^ \d,dD) . (4.6) 

Please observe, that the corresponding symbol for the trace space of i?°(d, D) will be 
2 .0(d, 91?) (and not VF2''^(d,Z3) as readers accustomed to the conventions used 

with Sobolev spaces might expect). 

In light of (|2.18p . the symplectic pairing on can be expressed as 

[l^,v]9d-= f ^^V, Lo,rjeW-^^\d,dD) . (4.7) 

JdD 

Whenever, W~^'^{d,dD) is treated as a real symplectic space, the pairing (14. 7p is 
assumed. The most important observation about (|4.7p is that the pairing [•, •] is utterly 
metric-free! 

Now we can specialize Theorem 13.81 to the curl operator. We give two equivalent 
versions, one for Euclidean vector proxies, the second for 1-fornis: 

Theorem 4.2 (GKN-theorem for curl, vector proxy version) . The mapping which 
associates to L C H^'^ (curlg, dD) the curl operator with domain 

©(curli) {v e H{curl,D) : 7t(v) e L} 

is a bijection between the set of complete Lagrangian subspaces of IjT^ (cvltIq^ dD) and 
the self-adjoint extensions o/curlmin. 

Theorem 4.3 (GKN-theorem for curl, version for 1-forms). The mapping which 
associates to L C W~'2'^(^d,dD) the -kd operator with domain 

V{-kdL) {v e W^(d,D) : i*r] e L} 

is a bijection between the set of complete Lagrangian subspaces of W^'^'^{d,dD) and 
the self-adjoint extensions of-kd defined on W^{d, D) . 

We point out that the constraint i*rj € L on traces amounts to imposing linear 
boundary conditions. In other words, the above theorems tell us, that self-adjoint ex- 
tensions of curlmin will be characterized by demanding particular boundary conditions 
for their argument vector fields, cf. [34]. 

Remark 4- Thanks to (j2.6p the symplectic pairing on W~^'^{d,dD) commutes 
with the puUback. Thus, if Z?, _D C M'^ are connected by a Lipschitz honiomorphism 
$ : I— > Z?, we find that $* : A^{dD) i-^ A^{dD) provides a bijective mapping 
between the complete Lagrangian subspaces of W~^'^{d,dD) and of W~~^'^{d,dD). 
Thus, pulling back the domain of a self-adjoint extension of curlmin on _D to D will 
give a valid domain for a self-adjoint extension of curlmin on D. In short, self-adjoint 
extensions of curlmin are invariant under bijective continuous transformations. This 
is a very special feature of curl, not shared, for instance, by the Laplacian —A. 

5. Hodge theory and consequences. Now we study particular subspaces of 
the trace space W~^'^{d, dD). We will take for granted a metric on dD that induces 
a Hodge operator 
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5.1. The Hodge decomposition. Let us first recall the well-known Hodge 
decomposition of spaces of square-integrable differential 1-forms on dD. For a more 
general exposition we refer to [29] : 

Lemma 5.1. We have the following decomposition, which is orthogonal w.r.t. the 
inner product of L^{h}{dD)): 

L^{A\dD) = dW°{d,dD)®irdW°{d,dD)®H\dD) . 

Here, T-0^{dD) designates the finite-dimensional space of harmonic 1-forms on dD: 

H^{dD) := {uj e L^{A^{dD)) : d w = and d *w = 0} . (5.1) 

In terms of Euclidean vector proxies, the space L'^{A^{dD) is modelled by the 
space Til{dD) of square integrable tangential vector fields on dD. Then, the decom- 
position of Lemma 15.11 reads 

mOD) = gradg H\dD) © cuTlgH^idD) © U^{dD) , 
V}{dD) := {v e l4{dD) : curlgv = and divgv = 0} . 

The Hodge decomposition can be extended to W^^'^{A,dD) on Lipschitz domains, 
as was demonstrated in [9, Sect. 5] and [5]. There the authors showed that, with a 
suitable extension of the surface differential operators, that 

H-5 (curia, az?) =gYa,Ai)H^{dD)®c\iY\9Hi{dD)®H^{dD) , (5.2) 

where, formally, 

Hi{dD) := Ag^H^^dD) , H:^^{dD) {v e H-i{dD) : [ vdS = 0} . 

JdD, 

(5.3) 

with dDi the connected components of dD. 

For C^-boundaries this space agrees with the trace space of H^{D). 

The result l|5.2p can be rephrased in the calculus of differential forms: 
Theorem 5.2 (Hodge decomposition of trace space). We have the following 

orthogonal decomposition 

W-i'\d,dD) = dW-^^°{d,dD)®*dwi^°{dD)®H^{dD) , (5.4) 

whith 

W^'°idD) Ag^{^ e W-^^°{dD : 1)90, - 0} . (5.5) 

with dDi the connected components of dD. 

The first subspace in the decomposition of Theorem 15.21 comprises only closed 
1-forms, because 

d(dW^-iO(d,9D)) =0. (5.6) 

The second subspace contains only co-closed 1-forms, since 

d*(*dM^i'°) =0. (5.7) 

The Hodge decomposition hinges on the choice of the Hodge operator 7k-. Conse- 
quently, it depends on the underlying metric on dD. 
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5.2. Lagrangian properties of the Hodge decomposition. We find that 
the subspaces contributing to the Hodge decomposition of Theorem 15.21 can be used 
a building blocks for (complete) Lagrangian subspaces of W~^'^{d, dD). 

Proposition 5.3. The linear space is a Lagrangian subspace of 

W~'^'^{d,dD) (w.r.t. the symplectic pairing -jgj^) 

Proof. We have to show that 

[dw,dry]g^ = Vw, ry e VF"'^"(d, 9D) . (5.8) 

By density, we need merely consider uj, rj in ^^"(d, dD). In this case, it is immediate 
from Stokes' Theorem (dD has no boundary) 

[diLj,di]]gD = / diuAdrj = / cj A d^ ry = . (5.9) 

JoD JdD 

□ 

Proposition 5.4. The linear space is a Lagrangian subspace of 

W-i^'^id, dD). 

Proof. The proof is the same as above, except that one has to use that * is an 
isometry with respect to the inner product induced by it: 

[kduj,*dri]dD = / *duj A^drj = - / duj Adrj = - / A d^ ?/ = . (5.10) 

JdD JdD JdD 

□ 

In a similar way we prove the next proposition. 

Proposition 5.5. The space of harmonic 1-forms Ti}{dD) is symplectically or- 
thogonal to dW-^'^{d,dD) and*dW^'°{dD). 

The Hodge decomposition of Theorem 15.21 offers a tool for the evaluation of the 
symplectic pairing [•, ■]gj^. Below, tag the three components of the Hodge decomposi- 
tion of Theorem 15.21 by subscripts 0, _L and H: for u!,ri ^ W~^'-^{d, dD) 

UJ = dujQ + -*dLu± + ujh and rj = d 770 + *dr]± + rjn ■ (5-11) 

Note that the forms wq and u)^_ are not unique since the kernels of d and 7k- d are not 
empty (they contain the piecewise constants on conected components of dD). 

Taking into account the symplectic orthogonalities stated in Propositions [531 [531 
and l5.5l 

f [dcjo,d77o]aD = [*duj^,-kdrn_]oD = [dujo,riH]dD = [★dcJi,?;_ff]aD 
\ = [u;H,di]c]dD = [uJH,-*:d'i]±]oD = 

(5.12) 

we see that we can compute the symplectic pairing on W~'^'^{d, dD) according to 

[uj,r]]9D = [dujQ,^dri±]9D + [*dujj_,dr]Q]QD + [uJH,r]H]dD ■ (5.13) 

It can also been expressed in terms of the L^-inner product (more precisely, its ex- 
tension to duality pairing) as 

[uj,ri]dD = idujQ,dr/±)i,9D - {duj±,dr]o)i^oD + [(^H,r]H]dD ■ (5-14) 
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5.3. The symplectic space H^{dD). Let us recall that the space of harmonic 
1-forms on dD (a 2 dimensional compact C°°-manifold without boundary) is a finite 
dimensional linear space with 

dim{n^{dD)) = 2g, (5.15) 

with g the genus of the boundary, that is, the first Betti number of D. The reader 
can refer to Theorem 5.1, Proposition 5.3.1 of [4] and Theorem 7.4.3 of [29]. 

Since the set of harmonic vector fields is stable with respect to the Hodge operator 
(note that = — 1 for 1-forms on dD) 

rieH\dD) =^ 7] e L^{A\dD)), d r/ = 0, d*r/ = 

^ *r] e L'^{A^idD)), d*(*ry) = 0, d(*?7) (5.16) 

Thus we find that the pairing [•, -J^^ is non-degenerate on TC^{dD): 

([i^H,riH]dD = 0, V?7_ff e H^{D)^ [uJH,-*^t^H]dD = {'^H-,t^H)l,dD = . (5.17) 

Lemma 5.6. The space of harmonic 1-forms Tl^{dD) is a symplectic space with 
finite dimension when equipped with the symplectic pairing [■,-]dD- It is a finite- 
dimensional symplectic subspace ofW^"S'^{d,dD). 

6. Some examples of self-adjoint curl operators. Starting from the Hodge 
decomposition of Theorem 15.21 we now identify important classes of self-adjoint ex- 
tensions of curl. We rely on a generic Riemannian metric on dD and the associated 
Hodge operator. 

6.1. Self-adjoint curl associated with closed traces. In this section we aim 
to characterize the complete Lagrangian subspaces L of 2'i(d, 9Z3) (equipped with 
[■j "la-o) which contain only closed forms: 

L C Z-i'\dD) -.^ {rj eW~i'\d,dD) : d?7 = 0}. (6.1) 

Hodge theory (see Theorem 15. 2p provides the tools to study these Lagrangian sub- 
spaces, since we have the following result: 

Lemma 6.1. The set of closed 1-forms in W~^'^{d,dD) admits the following 
direct (orthogonal) decomposition 

Z~-^'^{dD) = dW-i'"id,dD) ® li}{dD) . (6.2) 

Proof For uj e Z^i-^{dD), -kduj±_ of ([5lT|) satisfies 

d(*dw_L) = 0, d*(*dtj_L) = 0, {irdLj±)H = (6.3) 

which implies that ★ d lu± — and yields the assertion of the lemma. □ 

The next result is important, as it states a necessary condition for the existence 

of Lagrangian subspaces included in Z~2'^(dD). 

Lemma 6.2. The space Z~'S'^{dD) includes its symplectic orthogonal 

dW-^^°{d, dD). 
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Proof. Recall from Definition 13.21 that the symplectic orthogonal of Z^i'^{dD) is 
defined as the set 

{cu <EW-i^'^{d,dD) : [uj,f]]9D =0, ^/t] e Z-^'\dD)} . (6.4) 

Using Theorem 15.21 for u; = d + * d uj± + uir and Lemma 16.11 for rj — d r/Q + tjh, we 
have with (|5.13|) : 

[uj,ri]9D = [*rduj±,dT]a]gD + [uJH^VHho- (6.5) 

This impHes with 77 = -kuJH = Vh (here we use the stability of Ti}{dD) with respect 
to the Hodge operator) 

[u},-*^UJH]dD = [(^H,-*'LOH]dD = / ^if A -kUJH = =^ LUH ^0 , (6.6) 

JdD 

and, for 7] — d rjo with rja ^ uj^ & W'^/^'°{dD) 

[uj,duj±]Q£i = [kduj±,dLOj_]£)D = dLu±A*dLU±_ dLu±—0. (6.7) 

JdD 

Hence, we have oj — dcuo (and ujh = 0). The converse holds due to (|6.5p . □ 

Lemma 16.21 tells us that, when restricted to the subspace of closed forms, the 
bilinear pairing [-j-Jod becomes degenerate. More precisely, on the subset of closed 
forms, one can use the splitting (|6.5p and evaluate [■, ■]dD on Z~2'^(9D) according to 

[u;,'q]QD = [i^H,VH]dD , Vw, 77 e Z"3.i(ai:)) . (6.8) 

Hence, this pairing depends only on the harmonic components. Thus another mes- 
sage of Lemma [6.21 is that [•, ■]od furnishes a well-defined non-degenerate symplectic 
pairing, when considered on the co-homology factor space 

m\dD,R) = Z-^2'\dD))/ dW-^^°{d,dD) . (6.9) 

This space is algebraically, topologically and symplectically isomorphic to 7i^{dD), 
the space of harmonic 1-forms, see Section 15.31 

This means that all the complete Lagrangian subspaces L of W^^'^{d,dD) con- 
tained in 2^i(9£)) are related to complete Lagrangian subspaces L-h of H^{dD) (or 
equivalently to the complete Lagrangian subspace Lh of M^{dD,R)) by 

L = dW-^'"{d,dD) O Ln (or, equivalently, = L/ d T4^-^^°(d, 9D)) . 

(6.10) 

Thus, we have proved the following lemma (the symplectic pairing [^-Ja/j is used 
throughout) 

Lemma 6.3. There is a one-to-one correspondance between the complete La- 
grangian subspaces L of the symplectic space W~'^'^{d,dD) satisfying 

LcZ-^^\dD) (6.11) 

and the complete Lagrangian subspaces Lu of Ti}{dD). The bijection is given by 

Theorem 14.31 and Lemma [^751 lead to the characterization of the self-adjoint curl 
operators whose domains contain only functions with closed traces. 
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Theorem 6.4. There is a one-to-one corresondance between the set of all selfad- 
joint curl operators * ds satisfying 

V{* ds) C \^LU e W\d, D) i*uj e Z-i^\dD)^ (6.12) 

and the set of complete Lagrangian suhspaces L-j-c ofT-L^{dD). They are related accord- 
ing to 

V(*ds) = \^LueW\d,D) i*Lue dW-^'°{d,dD)®LHj ■ (6.13) 
Obviously, the constraint 

i*uj £ d W-^'°{d,dD) © Lh (6.14) 

is a boundary condition, since it involves only the boundary of the domain D. In 
addition, we point out that no metric concepts enter in ()6.13p . cf. Section IH 

Remark 5. Now, assume the domain D to feature trivial topology, that is, the 
genus of D is zero, and the space of harmonic forms is trivial. Theorem 16.41 reveals 
that there is only one self-adjoint -kd with domain containing only forms with closed 
traces 

V{*ds) = [lo eW\d,D) d(rw) = o|. (6.15) 

In terms of vector proxies, this leads to the self-adjoint curl operator with domain 

P(curl5) = |u e iy°(curl, D) curl(u) • n = on dD^ , (6.16) 

which has been investigated in [35,39]. In case D has non-trivial topology, then 
dim{H^{dD)) — 2g ^ and one has to examine the complete Lagrangian subspaces 
of H^{dD), which is postponed to Section [ 



6.2. Self-adjoint curl based on co-closed traces. In this section we seek 
to characterize those Lagrangian subspaces L of W^^'^{d, dD) that contain only co- 
closed forms, ie. 

L C |tJ e VK"^/^'^(d,9D) : d*tj = o|. (6.17) 

The developments are parallel to those of the previous section, because, as is illustrated 
by (|5.13p . from a symplectic point of view, the subspaces of closed and co-closed 1- 
forms occuring in the Hodge decomposition of Theorem 15.21 are symmetric. For the 
sake of completeness, we give the details, nevertheless. 

Lemma 6.5. The subspace of co- closed 1-forms of W^^^^'^{d,dD) admits the 
following orthogonal decomposition 

^^cu &W-^/^'\d,dD) : d*a) = o} = ir dW^/^'"{dD) © H\dD) . (6.18) 

Proof. For oj co-closed, d in (|5.1ip satisfies 

d(d ujo) = 0, d ★(d ujo) ^ 0, (d ljo)h =0 (6.19) 
which implies that dujo — and proves (|6.18p . □ 
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The next result is important as it states a necessary condition for the existence 
of Lagrangian subspaces comprising only co-closed forms. 

Lemma 6.6. The symplectic orthogonal of the subspace of co- closed forms of 
W-^/^^\d,dD) isi^dW^/^^^idD). 

Proof. Use the definition of the symplectic orthogonal as 

{uj e W-^/'^'\d,dD) : [uj,r]]9D = V?7 co-closed} . (6.20) 

Using Theorem 15.21 for lu — d ojq + d uj± + ujh and Lemma 16.51 for rj — *dr]± + rjn, 
gives 

['^,v]dD = [dt^o,*d77^]aD + [t^_ff,?7ff]a_D . (6.21) 

Choosing rj = -kujH = r/H (here we use the stability of H}{dD) with respect to the 
Hodge operator) this implies 

[u},-*:UJH]dD = [uJH,-*'LOH]dD = / ^if A -kUJH = LOh = 0, 

JdD 

and, for 77 = ★d 77^ with 77^ £ wi '^{dD) 
[uj,-^dr]±]dD = [dujQ,-kdrn_]sD = / dwoA*d7/_L = =^ d*(dwo)==0. 

JdD 

Moreover, one has d(d ujq) = and (d ujo)h = 0, which shows that divo — 0- 

Hence, we have w = *d (d wq = and uh = 0). The other inclusion holds due 
to fOTjl . □ 

Remark 6. Formally, in the proof we have used 77 = * d r/^ with ri±^ = lO{) which 
shows that 

[ui,r]\aD = [dcjo,*dwQ]a_D = / dwoA*dcjo =^ dcjo = 0. (6.22) 

JdD 

However, the lack of regularity of wq does not allow this straightforward computation. 

When restricted to the space of co-closed forms, the bilinear pairing [•,-]a_D be- 
comes degenerate. However, due to Lemma 16. 6[ it is a non-degenerate symplectic 
product on the co-homology factor space. 

^uo eW~^'^^^{d,dD): d*tj = 0}/*dW^3/2.0(^£,) ^ (g 23) 

which can be identified with H^{dD). Indeed, also on the subset of co-closed forms, 
one can evaluate [•,-]dD by means of (|6.8p . 

Lemma 6.7. The complete Lagrangian subspaces L ofW~^^'^'^{d,dD) containing 
only co-closed forms are one-to-one related to the complete Lagrangian subspaces Lq-i 
ofn\dD) by 

L = -kdW^/^^"{dD)®Ln ■ (6.24) 

Theorem 14.31 and Lemma 16.71 lead to the characterization of the self-adjoint curl 
operators based on coclosed forms: 

Theorem 6.8. There is a one to one correspondance between the set of all self- 
adjoint operators ★ ds satisfying 

V{-kds) c c,W^{d,dD) : di^{i*Lu) = Q^ (6.25) 
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and the set of complete Lagrangian subspaces Ln of Ti.^{dD) equipped with [-j-jdD- 
The underlying bijection is 

V{-kds) = [uj eW^{d,D) : i*Lje irdW^^^'°(dD)®Ln} ■ (6.26) 

Remark 7. Let 13 be a domain with trivial topology. Then there is only one self-adjoint 
operator * d whith domain containing only forms whose traces are coclosed 

V{*ds) = |w e iy"(d,17) : d*(i*w) = o|. (6.27) 

In terms of Euclidean vector proxies, we obtain the self-adjoint curl operator with 
domain 

P(curls) = ^^ue H{curl,D) : diva(7t(u)) = on (6.28) 

On the contrary, if D has non trivial topology, then one has to identify the complete 
Lagrangian subspaces of Ti}{dD). This is the topic of the next section. 

6.3. Complete Lagrangian subspaces of 7Y^(9D). The goal is to give a rather 
concrete description of the boundary conditions implied by (|6.13p and (|6.26p . Concepts 
from topolgy will be pivotal. 

To begin with we exploit a consequence of the long Mayer- Vietoris exact sequence 
in CO- homology [4] , namely the algebraic isomorphisms [23] 

n^{dD) '^m^{dD-M.) = i*„tf(D;]R) +CtEI^(^';K) . (6.29) 

Here, m^{D) is the co-homology space Z^{d,D)l AW^{A,D), and i-,^ : dD ^ D, 
iout : dD (-> D' stand for the canonical inclusion maps. We also point out that [23] 

^ dimM^ {dD;R) dimHi(£);M) = dimHi(I3';R) = g , (6.30) 

where g £ No is the genus of D. 

Next, we find bases oiM^{D; R) and M^{D'; MVusing the Poincare duality between 
co-homology spaces and relative homology spaceqj 



M^{D- R) ^ M2{D, dD; R) . (6.31) 
Consider the relative homology groups (with coefficients in Z) 

m2iD,dD;Z) and m2{D\dD;Z) (6.32) 
as integer lattices in the vector spaces 

m2{D,dD;m) and W2iD',dD;R). (6.33) 



In this article, we denote by 

• BIi(j4; R) the i*'' homology group of A with coefficients in R; 

• B1*(j4; R) the i*'' co-homology space of A with coefficients in R; 

• Hi (A, B; R) the i"" relative homology group of A relative to B with coefficients in _R; 

• H' (j4, B; R) the i"" relative co-homology space of A relative to B with coefficients in R. 
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In [21], it is shown that these lattices as Abehan groups are torsion free, and that 
representatives of homology classes can be realized as orientable embedded surfaces. 
More precisely, we can find 2g compact orientable embedded Seifert surfaces, or "cuts" 

S,,Sl, Is^is^s, (6.34) 

such that their equivalence classes under appropriate homology relations form bases 
for the following associated lattice^ 

{(^^>}Li for Hap, a/?; Z), {{S'^}}'^^, for M^iD' , dD;Z) (6.35) 

and 

{{dSl),{dS^)}l^, forHi(9D;Z). (6.36) 

In other words, the boundaries dSi, dS[ provide fundamental non-bounding cycles on 
dD. 

In [23], it was established that the set of surfaces {(5'i)}f^j U {(S-)}^^! can be 
chosen so that they are "dual to each other" . Here, this duality is expressed through 
the intersection numbers of their boundaries, see Chapter 5 of [18]. 

Lemma 6.9. The set of surfaces {{Si)}l^i U {{S'^Yi=i ca"' &e chosen such that 
the intersection pairing on Wi{dD]'L) can he reduced to (I <i,j < g) 

Int((9S,),(55j)) = <5,,, , 

Int((9S:),(a5,)) = -<5.,,. ^''■■"> 

Furthermore, when the boundaries of these surfaces are "pushed out" of their 
respective regions of definition, wc get curves in the complementary region 

dSl — > a, ^S^ — > C[. (6.38) 

The homology classes of these curves form bases for homology lattices as follows 

{{C^)]Ll forHi(Z?;Z), and {(C^ILi for Hi(i?'; Z). (6.39) 

This paves the way for a construction of bases of the co-homology spaces on D and 
D' [23]: 

Lemma 6.10. The co-homology classes generated by the closed 1-form in the sets 
defined for 1 < i < g 

e L'^{A^{D)) : dni = and J m ^ Sij for 1 < j < g| 

e L^{A\D' j) : dn^^O and j n[ ^ (5y for 1 < j < 5} 

form bases ofM^{D;M.) and M.^ {D' ;M.) , respectively. 

For instance, Hi can be obtained as the piecewise exterior derivative of a 0-form 
on £> \ Si that has a jump of height 1 across Si. An analoguous statement holds for 
k'^ with Si replaced with 5*^'. More precisely, one has for 1 < i < g 

3^, £ VF°(d, D) : = d V'j on D \ 5 and [V-^ls, = S^j (6.40) 



^Throughout the paper (•) denotes the operation of taking the (relative) homology class of a 
cycle 
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3Vj^ e W°{d,D) : = d V- on D' \ S' and [i/;^]s^ = S,^, (6.41) 

with [-Jr denoting the jump across T. 

Lemma 6.11. For 1 ^ m,n ^ g, we have 

a) / i*„{Krn) M*„{Kn) = 0, 
JdD 

(6.42) 

b) I Ct(«™)A«Lt«) = 0. 

Proof. To estabhsh a) we rewrite the integral as one over D, as the foUowing 
calculation shows 



dD JdD JD 

(d Krn) A K„ - A (d K„) = 0. 



D 

Similarly, b) follows from an analogous calculation where dD = —dD' with forms 
defined on _D'. □ 

Lemma 6.12. For 1 ^ i, j < g, we have 

Kni^i) f^^uti^j) = kr (6.43) 

dD 

Proof. Let us represent the 1-form m by means the 0-form ipi, which jumps across 
Si, see (|6.40p . Taking into account that d il^t'^'j — 0> S^^ 

^^nl^^ ^ i*out^'j = cl i-^.V'j A = d (i*„?A, A . (6.44) 

Applying Stokes Theorem leads to (one has to take care of the orientation) 

im'«»ACt'«j= / [CV'AC^k'] + / Ml^tHi'] . (6.45) 

dD JdSi ' JOS', 



By ([OO)) . we get 



A = / + 0. (6.46) 

9D JdSi 



Since 95^ G (C'j'); the result follows from (|6.38p and Lemma [61101 □ 

Remark 8. When the cuts do not satisfy (|6.37p , a generalization of Lemma 16.121 
takes the form 

A »:„*(«;■)- Int((a5,),(9^j)). (6.47) 

dD 

with \ni[{dSi), (dSj )) the intersection number of {dSi) and (dSj), see [18]. 

Now, take (|637|) for granted. Write f^H.i, f^'n ij 1 * ^ 5i for the unique harmonic 
1-forms, i.e., nn^iTi^'fji 

e n^{dD), such that 
i*^Ki = KH,i + d a , i*ut«^i = i^'h,! + d /? , (6.48) 
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for some a, (3 ^ L'^{hP{dD)). Combining Lemmas 16.111 and 16.121 gives the desired 
symplectic basis of the space of harmonic 1-forms on dD: 

Lemma 6.13. The set {i^H,iT >^'h i}l=i symplectic basis ofH}{dD). 

Obviously, since the trace preserves integrals and integrating a closed form over 
a cycle evaluates to zero, the 1-forms kra and n'jj ^ inherit the integral values over 
fundamental cycles from Ki and k^, c/. Lemma 16.101 

HH,i^5tj, I KH,i^O, / K'jjt^Sij, / K^j = 0. (6.49) 
aSj Jas'. Jdsr ' JdSj 

Lemma 6.14. Given interior and exterior Seifert surfaces 5*^,5'^', the condi- 
tions (j6.49p uniquely determine a symplectic basis {kh,i, ■ • • i >^H,gj i^'h ii ■ ■ • j '^'h g) 
n^{dD). 

Proof. If there was another basis complying with (|6.49p . the differences of the 
basis forms would harmonic 1-forms with vanishing integral over any cycle. They 
must vanish identically. □ 

Given a symplectic basis, we can embark on the canonical construction of complete 
Lagrangian subspaces of H^{dD) presented in Remark [31 We start from a partition 

lUl' = {l,...,g} , /n/' = 0. (6.50) 

Owing to Lemma 16.121 and ()6.37|) the symplectic pairing [•, ■]gD has the matrix repre- 
sentation 



Ogxg Igxg 
-Igxg Ogxg^ 



2g 



(6.51) 



with respect to the basis 

({tff.ilie/ U {-k'h, JjG/') U ({-k'^, JjG/ U {KH,i}iei'^ (6.52) 
oin^idD). Thus, 

L-H ■■= span{KH,t}iei U {-K^^ijiei' (6.53) 

will yield a complete Lagrangian subspace of H^{dD). By theorems 16.41 and 16. 8[ L-j-c 
induces self-adjoint curl — *d operators. From Lemma Lemma [^31 and (|6.49p we 
learn that their domains allow the characterization 



V{curls) := {uj eW\d,D) : d{i*u;) = 0, [ uj = 0,jel,[ u; = 0,jel'} 

JaSj Jas'^ 

(6.54) 



in the case of closed traces, and 



I?(curl,) := {w G l^^(d,D) : d*(i*cj) =0, / w = 0, j G /, / = 0, j G /'} , 

JaSj Jas'^ 

(6.55) 

in the case of co-closed traces, respectively. In fact, the choice /' = together with 
closed trace is the one proposed in [39] to obtain a self-adjoint curL 
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7. Spectral properties. Having constructed self-adjoint versions of the curl 
operator, we go on to verify whether their essential spectrum is confined to and 
their eigenfunctions can form a complete orthonormal system in L^{D). These are 
common important features of self-adjoint partial differential operators. 

The following compact embedding result is instrumental in investigating the spec- 
trum of curls. Related results can be found in [37] and [32]. 

Theorem 7.1 (Compact embedding). The spaces, endowed with the W^{d,D)- 
norm, 

Xo :={lu e W\d, D) : d* uj = 0, i*{*uj)^0} 
and :={w eW^{d,D) : d" uj = 0, di.{i*uj) = 0} 

are compactly embedded into L'^{A^{D)). 

Remark 9. In terms of Euclidean vector proxies these spaces read 

Xo ={v e if (curl, D) : divv = 0, 7„u = 0} , 
X-^ ={v e if (curl,!)) : divv = 0, diva(7fu) = 0} 

where the constraint diva(7fu) = should be read as "orthogonality" to 
gradg H^{dD) in the sense of the Hodge decomposition. 

Proof, [of Thm. 17. 1| The proof will be given for X-^ only. The simpler case of Xq 
draws on the same ideas. We are using vector proxy notation, because the proof takes 
us beyond the calculus of differential forms. Note that the inner product chosen for 
the vector proxies does not affect the statement of the theorem. 

A key tool is the so-called regular decomposition theorem that was discovered 
in [3], consult [19, Sect. 2.4] for a comprehensive presentation including proofs. It 
asserts that there is C > depending only on D such that for all u e J? (curl, D) 
there are functions # £ {H^{D))^, (p e H^{D), with 

ff(curl,_D) • 

(7.1) 

Let (u„)jjgpj be a bounded sequence in X^ that is 

divu„ =0 in D and diva(7tu„) = on dD , (7.2) 
3O0: ||u„||^.(^) + ||curlu„||^.(^) <C. (7.3) 

Write u„ = + grad(/p„ for the regular decomposition according to (|7.ip . Thus, 
($„)„gN is bounded in {H^{D))^ and, by ReUich's theorem, will possess a sub- 
sequence that converges in L^{D). We pick the corresponding sub-sequence of (u„)„gN 
without changing the notation. 
Further, 

divu„ = ^ -A(^„ = div*„ (bounded in £^(1?)) , (7.4) 
diva(7tu) = =^ -AQoh^Pn) ^ divailt^n) (bounded in i/"^ (91))) . (7.5) 

We conclude that {j(pn)neN is bounded in H^{dD) and, hence, has a convergente 
sub-sequence in H2{dD) (for which we still use the same notation). The harmonic 
extensions (pn of jipn will converge in H^(D). 

Finally, the solutions ^„ G H^{D) of the boundary value problems 

-A(^„=divf>„ in D , (pn = on dD , (7.6) 
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will possess a sub-sequence that converges in H^{D), as (— Aoir) ^L'^{D) is compactly 
embedded in H^{D). Since ipn = ^„ + ^r! , this provides convergence of a subsequence 
of (*„ + grad^„)„gpj in L^{D). □ 

Let curls : T^s C L'^{A^{D)) L^(A^(_D)) be one of the self-adjoint realizations 
of curl discussed in the previous section. Recall that we pursued two constructions 
based on closed and co-closed traces, respectively. 

Remark 10. Even if the domain 2?s of the self-adjoint curls is known only up to 
the contribution of a Lagrangian subspace of we can already single out special 
subspaces of 7)^: 

(i) For the curl operators based on closed traces, see Sect. 16.11 in particular 
Thm. 16.41 we find 

AW°{d,D)(zVs. (7.7) 

Indeed, for G d VF°(d, D) there exists -q G Vl^"(d, D) with uj ~ drj. Due to 
the trace theorem, i*ri belongs to W'^^ {A,dD). Consequently, it follows from 
the commutative relation (|2.6p that i*uo = di*ri belongs to dW^^ {d,dD). 
We conclude using (|6.13p . 

(ii) For the curl operators based on co-closed traces introduced in Sect. 16.21 it 
follows that 

dW^id,D)cVs. (7.8) 

This is immediate from the fact that 

r] G Wg (d, £)) and uj = dr] implies i*oj = d i*r] — , (7.9) 

which means that uj belongs to Vg, see (I6.26p . 
In the sequel, the kernel of curls will be required. We recall that 

A/'(curls) = V, n A/'(curlmax) 

is a closed subspace of _L^(A^(_D)). Moreover, since d^ = and due to (|7.7|) and (|7.8|) . 
one has 

d W"{d, D) C 7V(curls) in the closed case, (7.10) 
d W^{d.D) C A/'(curls) in the co-closed case. (7-11) 

Lemma 7.2. The operator curls is hounded from below on Vg n A/'(curls)^; 

3C = C{D): < C||curlsw|| Vo) G Ps n A/'(curls)^ . 

Proof. The indirect proof will be elaborated for the case of co-closed traces only. 
The same approach will work for closed traces. 

We assume that there is a sequence (ijj„)^j^gj^ C nA/'(curls)^ such that 

||w„|| = l , ||curlw„|| < VnGN. (7.12) 

Since (jJ„ G A/'(curls)^, the inclusion (|7.1ip implies that = 0. As a consequence 
of (|7.12p . (i^n)„gN is a bounded sequence in X^. Theorem 17.11 tells us that it will 
possess a subsequence that converges in _L^(A^(_D)), again we call it (aj„)^gpj. Thanks 
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to (|7.12p it will converge in the graph norm on Vg and the non-zero limit will belong 
to 7V(curls) n7V(curls)^ = {0}. This contradicts ||tj„|| = 1. □ 

From Lemma [7^ wc conclude that the range space 7?.(curls) is a closed subspace 
of i^(A^(_D)), which means, 

7^(curl^) = 7V(curl,)^ . (7.13) 

Thus, we are lead to consider the symmetric, bijective operator 

C := curl^ : V, n Af{curh)^ C Af{curh)^ k-> 7V(curl,)^ . (7.14) 

It is an isomorphism, when T>s H A/'(curls)^ is equipped with the graph norm, and 
A/'(curls)-'- with the L^(A^(Z?))-norm. Its inverse is a bounded, self-adjoint op- 
erator. 

Theorem 7.3. The operator curl., has a pure point spectrum with oo as sole 
accumulation point. It possesses a complete -orthonormal system of eigenfunctions. 
Proof. The inverse operator 

Q-^ : AA(curl,)^ ^V^n AA(curl,)^ (7.15) 

is even compact as a mapping i^(A^(_D)) L'^{A^{D)). Indeed, due to (|7.10p and 
(|7.1ip the range of C^^ satisfies 

Vg nA/'(curls)^ C Xa in the closed case, (7-16) 
Vs n7V(curls)-'- C in the co-closed case. (7.17) 

By Theorem 1 7. 1[ the compactness follows. Riesz-Schauder theory [40, Sect. X.5] tells 
us that, except for its spectrum will be a pure (discrete) point spectrum with 
zero as accumulation point and it will possess a complete orthonormal system of 
eigenfunctions. 

The formula, see [38, Thm. 5.10], 

X-^ - Q-^ = X-\C - X)C-^ (7.18) 

shows that for A 7^ 0, 

• A^^ — C^^ bijective => C — A bijective , 

• A/'(A-i -C-i) =AA(C-A) . 

Thus, cr(C) = {a-{C^^) \ {0})^^ and the eigenfunctions are the same. □ 
8. curl and curl curl. 

8.1. Self-adjoint curl curl operators. In the context of electromagnetism we 
mainly encounter the self-adjoint operator curl curl. Now we explore its relationship 
with the curl operators discussed before. A metric on D and an associated Hodge 
operator * will be taken for granted. 

Definition 8.1. A linear operator S : V{S) c L'^{A''-{D)) ^ L'^{A'-{D)) is a 
curl curl operator, if and only if S is a closed extension of the operator -k d -k d defined 
for smooth compactly supported 1-forms. 
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Two important extensions of the curl curl operator are the maximal and the 

minimal c;xtcnsions: 

Lemma 8.2. The domain of the minimal closed extension (curl cur l)iiiin of the 
curl curl operator is 

Vrain = {w G W^{d,D) : *dio e W^{d,D)^ (8.1) 

or, equivalently, in terms of Euclidean vector proxies 

Vmin =1 u e L^(D) : curlu e i^(D), curl curl u e L'^{D), 
7f (u) = 0, and 7t(curl (u)) = on dD^. 

The adjoint o/ (curl curl),„in is the maximal closed extension (curl curl) max- It is 
an extension of the curl curl operator with domain 

I'max = Vi®V2, (8.2) 

with 

={w e WqH £/,£>) : *duj€ W^{d,D)y (8.3) 

V2 =|w e L^{A\D)) : c/*c/w = o| . (8.4) 

Proof. The domain fj^jj^ of the minimal closure is straightforward. We recall the 
definition of the domain of the adjoint T* of an operator T : I'(T) C H 1-^ H 



V 



(T*) = {ueif : 3C„ >0: {u,Tv)h ^ \\v\\h V« G I?(T)} . (8.5) 



Let I^max stand for the domain of the adjoint of the minimal curl curl operator. First 
we show that 

ViQV^C Pma^ . (8.6) 

Let us consider lj € Vi and r] G "Dmin- By integration by parts and the isometry 
properties of * we get 

y^Ad^d?? = y d^do; A?7 < ||d*da;|| IItjII . (8.7) 

D D 

This involves T>i C X>max- 

Now we consider a; G 2?2- The relation d*da; = has to be understood as 

Jd-kdLjAri = V?7 G A^(£)) smooth, compactly supported . (8.8) 

D 

As the smooth compactly supported 1-forms are dense in Dmin with respect to the 
topology induced by the norm 

ILII + II curl(a;)|| + || curl(curl(a;))|| , (8.9) 



Self-adjoint curl operators 



25 



it follows that 

JujAdi.dri = Vr/ e P^in , (8.10) 

D 

and, finally, X'2 C Pmax- This confirms (|8.6p . 
Next, we prove 

2?max C Dl ©I?2 . (8.11) 
Pick, uj 6 I'max- There exists ip G i^(A^(£))) such that 



LuAd-kdri;= J ip A-kr] V77 £ Pmin • (8-12) 

D D 

Since d* (yS = (pick r] ^ diy in (I8.12p ). and j^^p A = for j-jn e H^{D), there 
exists wi e M^^(d,£)) satisfying 



★ d*da;i = (y3 in D, 
i*0Ji ~ on dD. 



.13) 



Note that this lo\ belongs to T>\. Then ^2 = — ^^i satisfies 

J{uj-uJi) Adirdri = Vr/ e P^in =^ d*dw2=0. (8.14) 

D 

It follows that UJ2 G I?2 • Since = tJi + (^2 , we have proven (|8.1ip . □ 
Remark 11. The last lemma gives a nice example for 

(T^)* ^ (T^)*. 

Indeed, the minimal extension of the formal curl curl boils down to the squared 
minimal curl operator curlmin with domain W^{d,D) 

(curlcurl)min = curlmin CUrlmin 

The adjoint of curlmin is the curlmax operator with domain T/F^(d, D), but 
(curl curl),„ax ^ curl,„ax curlmax • 



To identify self-adjoint curl curl operators we could also rely on the toolkit of 
symplectic algebra, using the metric-dependent symplectic pairing 

[uj,r]] = / d^dwAr?- / tjAd*dry. (8.15) 
Jd Jd 

As before, complete Lagrangian subspaces will give us self-adjoint extensions of 
(curlcurl)i„in that are restrictions of (curl curl)niax- However, we will not pursue 
this further. 
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There are two classical self-adjoint curl curl operators that play a central role in 
electromagnetic boundary value problems. Their domains are 

P((curlcurl)Dh) ={w G W^id^D) : *dw e W\d,D)y (8.16) 

P ((curl curl) Ncu) ={w e W\d,D) : *dw e W^{d,D)^ . (8.17) 

Both can be written as the product of a curl operator and its adjoint: 

(curl curl)Dir = curlinax curlmin , (curl CUrl)Nou = curlmin CUrl,nax • (8.18) 

Less familiar self-adjoint curl curl operators will emerge from taking the square of a 
self-adjoint curl operator as introduced in Section [B] 

8.2. Square roots of curl curl operators. It is natural to ask whether any 
self-adjoint curl curl operator can be obtained as the square of a self-adjoint curl. We 
start with reviewing the abstract theory of square roots of operators, see [38, Sect. 7.3]. 

Let S be a positive (unbounded) self-adjoint operator on the Hilbert space H. We 
recall from [38, Thm. 7.20] that there exists a unique self-adjoint positive (unbounded) 
operator R saytisfying 

S = R^ i.e. D(S) = D(R2) := {u e D(R) / Rue D(R)} and Su = R^u if u e D(S) . 

(8.19) 

Lemma 8.3 (domain of square roots). Let Ri and R2 be two closed densely defined 
unbounded operators on H with domains D(Ri), D(R2) C H. 
//R* Ri = R2 R2, that is, 

D{Rl Ri) = D(R; R2) and Vit G D(RJ Ri), R^ Rm = R^ R2U , 

then D(Ri) = D(R2). 

Proof. For i = 1,2, D{Ri) equipped with the scalar product (u,v)i = (u,v)h -I- 
(RiU, Riv)H is a Hilbert space. 

Let us first prove that D(R*Ri) is dense in D(Ri) with respect to (•, We consider 
u e D(R* R,)^ 

VvgD(R*RO, = (u,v), = (u,v)h + (R,u,R,v)h = (u,v + R* R,v)h (8.20) 

As Id -I- R* Ri is surjective from D(R* R^) to H, see [36, Theorem 13.31], u is equal to 
zero. 

Hence, the spaces D(Ri) and D(R2) share the dense subspace D(R* Ri) = 
D(R2 R2). Moreover, their scalar products coincide on this subset: 

(u,v)h-H(RiU,Riv)h = {u,v+ R* Riv)h = (u,v + R^ R2v)h = (u,v)h + (R2U,R2v)h. 

We conclude using Cauchy sequences. □ 

Surprisingly, the simple self-adjoint operator (curl curl)Dir does not have a square 
root that is a self-adjoint curl: 

Lemma 8.4. The curl curl operator curl^ax curlmin does not have a square root 
that is a self-adjoint curl. 

Proof. Let us suppose that T — curl^ax curlmin admits a curl self-adjoint square 
root S which implies that 

curlinax curl*^^^ = CUrlmax CUrlmin = T = = S S* . (8.21) 
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since curl^ax and curl„iin are adjoint and S is self-adjoint. Due to lemma [H31 we 
have D(curli„ax) = D(S) and therefore 

S = curl^ax (8.22) 

since S and curlmax are both curl operators. Clearly, this is not possible since curlmax 
is not self-adjoint. □ 

Remark 12. We remark that the same arguments apply to the operator 
(curlcurl)Ncu- 

8.3. curl curl 7^ curl curl* is possible. Finally, we would like to show that not 
all the self-adjoint curl curl operators are of the form R R* with R a curl operator. 

Following an idea of Everitt and Markus — a similar construction or the Laplacian 
is introduced in [16] — we consider the self-adjoint curl curl operator 

T" : P(T°) C -L^(L') I — >L^{D), ui — >curlcurlu (8.23) 

with domain 

V{T°) =. P„un ® 2?2 , (8.24) 

where Pmin and 2?2 are defined in (18. ip and (|8.4p . 

Proposition 8.5. There exists no curl operator R such that 

= R R* . (8.25) 



Proof. Suppose that there exists a curl operator R satisfying (|8.25p . By definition 
of the composition of operators one has 

^(T") = |u e D(R*) : Ru G D(R)| . 

Hence, this implies 

2?2CP(T") c D(R*) c W\d,D). 

This is not possible since 1^2 is not a subspace W^^(d, D). 

This can be illustrated by means of vector proxies and in the case of the unit 
sphere D. Consider the function 

+ OC 

u{r,9,z) = sinn6'^ , 

ri=l 

given the cylindrical coordinates. The curl and curl curl of u are 

+00 +00 
curlu = (^^^n'r"~^ cosnO^ er — (^"^^ n r'^~^ sin ri6'^ eg , 

n— 1 n— 1 

curl curl u = . 

Direct computation leads to 

llulP < +00 and ||curlu|| =+00. 



28 



R. Hiptmair, P.R. Kotiuga, S. Tordeux 



Hence, this u satisfies u £T>2 but u ^ Ji'(curl, D). □ 

Remark 13. In the same way, we show that there exists no curl operators Ri and 
R2 satisfying 

T° = Ri R2 . (8.26) 



Appendix. Frequently used notations: 

D bounded (open) Lipschitz domain in afRne space 

D' (compactified) complement D' := \ l) 

dD boundary of D 

n exterior unit normal vector field on dD 

u, V, . . . vector fields on a three-dimensional domain 

co,r],. . . differential forms 

V, u elements of a factor space/trace space of vector proxies 

Euclidean inner product in 

X cross product of vectors G 

T, S, . . . (unbounded) linear operators on a Hilbert space 

T* adjoint operator 

T^in The minimal closure of T 

Tmax The maximal closure of T 

25(7) domain of definition of the linear operator T 

Af{J) kernel (null space) of linear operator T 

TZ{T) range space of an operator T 

C°° (-D) space of infinite differentiable functions on D 

C°°{D) space of smooth vector fields {C°°{D)f 

Cq^{D) functions in C°°{D) with compact support in D 

C^Id) vector fields in {C^{D)f 

L^{D) real Hilbert space of square integrable functions on D 

L'^{D) square integrable vector fields in {L'^{D))^ 

i?(curl, D) real Hilbert space {v G L^{D) : curlv € L^{D)} with graph norm 

i?o(curl, D) closure of Cif(L») in if (curl, D) 

7t tangential boundary trace of a vector field 

7„ normal component trace of a vector field 

gradg surface gradient 

curia scalar valued surface rotation 

diva surface divergence 

d exterior derivative of differential forms 

A'^(M) differential fc- forms on manifold M 

A exterior product of differential forms 

■kg Hodge operator induced by metric g 

(•, ■)k,M inner product on A^{M) induced by a Hodge operator ★ 

I/^(A'^(M)) Hilbert space of square integrable fc- forms on M 

IHI normof L2(A'=(M)) ("L^-norm"): ||wf := {oj,uj)k,M 

W''{d,D) Sobolev space of square integrable fc-forms with square integrable 
exterior derivative 

W(f (d, D) completion of compactly supported A;-forms in W''{d, D) 

i* natural trace operator for differential forms 

[•, •] generic symplectic pairing 
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symplectic pairing of 1-forms on 2-manifold M 
jump of trace of a function across 2-manifold F 



L" symplectic orthogonal of subspace L of a symplectic space 

(•) (relative) homology class of a cycle 
Hi (A; R) homology group of A with coefficients in R 

]HI*(A; R) i*'* co-homology space of A with coefficients in R 

M.i{A, B; R) i*^ relative homology group of A relative to B with coefficients in 
R 

M^{A,B;R) relative co-homology space of A relative to B with coefficients in 
R 

'H}{dD) co-homology space of harmonic 1-forms on dD 

'M}{dD) first co-homology factor space of non-exact closed 1-forms on dD 

'Ri{dD) first homology factor space of non-bounding cycles on dD 

{■) selects (relative) homology class of a cycle 

W-^-^{6,dD) trace space of VKi(d,I?) 

Hj. ^ (curia, dD) tangential traces of vector fields in if (curl, D) 

Z~^{dD) closed 1-forms in W~^'^{d,dD) 

u^, u)-^ components of the Hodge decomposition of a; e W~^'^{d, dD) 
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